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. We study Wigner function value statistics of classically chaotic quantum maps on com- 

• pact 2D phase space. We show that the Wigner function statistics of a random state is 

^SJ ' a, Gaussian, with the mean value becoming negligible compared to the width in the semi- 

classical limit. Using numerical example of quantized sawtooth map we demonstrate that 
the relaxation of time-dependent Wigner function statistics, starting from a coherent initial 
state, takes place on a logarithmically short (cx log fi) time scale. 



The Wigner function^^ (WF) is an essential concept of phase space represen- 
tation of quantum mechanics, namely it is a useful and faithful representation of 
I a pure or mixed quantum state in terms of functions of canonical classical phase 

' space variables. It has many applications in various branches of physics, in par- 

ticular in quantum optics. However, WF cannot be interpreted as quantum phase 
I space distribution as it can develop negative values, in particular due to well known 

' oscillatory interference fringes following e.g. coherent wave-packet superpositions.'^^ 

WF played an important role in the realm of quantum chaology.^^'^-* In particular. 
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I it has been conjectured that WF of stationary eigenstates of bounded dynamical 

II ' systems in quasi-classical regimes condense, within classical scales, onto classically 

^ . invariant components of phase space. On the other hand, within a smaller - quan- 

ta I tum scale, e.g. of Planck cell size, the phase space structure of WF of typical, say 

random or ergodic states, is very much unknown. This question is important for 
the understanding of decoherence and quantum stability with respect to system's 
perturbations as discussed recently.^^'^)"'') 
^ ' In this paper we address the question of the structure of typical WF from a 

. statistical point of view. We define and analyze the statistical distribution of values 

of WF of a given quantum state. In particular, we are interested in the relaxation 
of this distribution with time, when we start from initial coherent state, and in 
the corresponding time scale. Here we limit ourselves on classically chaotic discrete 
time systems, namely chaotic quantum maps, where the full phase space is the only 
topologically transitive classical ergodic component. 

We show that for fully chaotic quantum maps the limiting WF value distribution 
is a Gaussian. The average value is fixed by normalization, while the second moment 
(or the variance) diverges in the quasi-classical limit (h 0), so we have roughly 
a symmetric distribution of positive and negative values of WF for random states. 
In addition, we show that for chaotic systems, the relaxation to equilibrium in a 
statistics of WF happens on a short oc log h (Ehrenfest) time scale. 
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We consider Wigner functions^-*' of a pure state from the Hilbert space 
corresponding to a classical system with a compact phase space x- For explicit 
calculation we choose 2D toroidal phase space x = "^2 that wc quantize following 
Rcf. 8). Let us write canonical Hilbert space basis of position states as |n), n = 
0, . . . , N — 1. Due to technical reasons is odd. WF is defined on a regular grid of 
N X N phase-space points 

W^{n,m) = ^^^e-2-"WiVj2^-2n+n'(^ + n'|V')(V'IO, 4 = ^1^6--'*^/^, 

^ n',l m' 

where all indices run from — (A^ — l)/2 to (A^ — l)/2, as will be assumed whenever 
we refer to the quantized torus. We study WF value distribution P{w) defined as 

P{w) = ^ j 5{w- W{x))dx, (1) 

where V is the phase space volume (area). In our case: (1/F) = ^. Let 

(. . .) denote the average with respect to P{w). The first two moments of P{w), the 
mean value W and the standard deviation a, are determined by the normalization 
of probability and purity = p of the state p = We shall 'normalize' WF 

by unitarizing the standard deviation 

cr^ = l and = (iV - 1)"^/^ (2) 

Particularly interesting is the relative standard deviation defined as k = a/W that 
scales with dimension N as k VN, i.e. it diverges in the semi-classical limit. 

The states of developed quantum chaos are associated with ergodic or random 
wave functions^) that can be defined in Hilbert space basis \l), I = 0, . . . ,N — 1 as 

N N 

1=1 1=1 

where q G C are random coefficients, that may be considered for large N as indepen- 
dent Gaussian variables^^^ of variance (|qP) = (. . .) denotes ensemble average 
over random states (3). WF value statistics of a random state at fixed point in 
phase space should be transitionally invariant and thus equal to the expected statis- 
tics of a 'chaotic Wigner function'. Therefore we can compute P{w) by analyzing 
the distribution of WF values at one phase space point, e.g 1^(0, 0), 



WiO, 0) = [ hi' + E + OcFq, I . (4) 



The expression (4) is a sum of N statistically independent terms, where the first 
I Cop is strictly non- negative and others have vanishing mean and finite fluctuation. 
The central limit theorem, for large N 1, implies that the distribution of W{0, 0) 
is a Gaussian 

^'ra.idom(^«) = exp Wf^ , (5) 
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with the first and the second moment equal to those of WF value distribution over 
the whole phase space (2). 

The dynamics is studied on the generic example of the quantized saw-tooth map 
on the toroidal phase space. The latter is with kicking strength K defined by the 
following evolution operator^^) 

[/ = exp I -i-p I exp ( i—n 

where n|n) = n\n), eyip{2iTip/N)\n) = |n + 1), and T = 2'k/N. Here we consider 
the regime > for which the classical sawtooth map is ergodic and uniformly 
hyperbolic. In the fig. 1 we show distribution P[w) of time-evolving state in the 
sawtooth map for K = 0.5 and A'^ = 2187, starting from a coherent state. We see that 
P{w) quickly relaxes from a particular initial form to the Gaussian distribution (5) 
predicted for random waves. The relaxation is the result of the wave self-interference 
caused by the corresponding classical stretching of the wave packet over the phase 
space. 
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Fig. 1. Time evolution of the WF value distribution (1) of the quantized sawtooth map for iiT = 0.5 
and N = 2187. Initial state is a coherent wave packet. Dashed curves give a theoretical Gaussian 
(5). 



To measure the deviation of distribution P{w) from the Gaussian statistics we 
define the excess e of the distribution by formula 

e = {w-WY/a"^ - 3. (6) 

Note, that e = for a Gaussian (5) and the size of e measures a deviation from Gaus- 
sian. Using the excess e we can quantitatively measure the transition of statistics to 
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the asymptotic Gaussian form, that happens on the Ehrenfest time scale 

tr ~ log(iV)/A, (7) 

where A is an effective Lyapunov exponent and N is the dimension of Hilbert space. 
Numerical experiments for the family of sawtooth maps show that tj- (7) has the 
right scaling, both with A*" (see fig. 2) and with A (see fig. 3). To check the rate 
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Fig. 2. The excess e of the WF value distrib- 
tion as a function of scaled time t/\ogN 
for several different values of N . The pa- 
rameter of the quantized sawtooth map is 
K = 0.5, and the initial state is a coherent 
wave packet. 
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Fig. ,3. Scaled time evolution of the excess e 
for the quantized saw-tooth map at dif- 
ferent kicking strengths K and fixed di- 
mension N = 2187. The Lyapunov ex- 
ponent is given by X{K) = log((2 + K + 
V(2 + i^)2-4)/2).") 



of breakdown of direct classical-quantum correspondence we propose to study the 
percentage of phase space with negative valued WF 

P-= [ P{w)dw (8) 

as a function of time. For random states P_ goes to 1/2 in the limit — > oo. In 
the fig. 4 we show the evolution of P_ in the quantized saw-tooth map, initially in a 
coherent state, for different dimensions N. We see that P_ converges very fast to 1/2 
on a time scale tc proportional to log(A'') but significantly smaller than relaxation 
time tr- However, both time scales have the same scaling with N and A. 

We have studied WF statistics of random states, and of time-dependent states 
of time-periodic, e.g. kicked quantum systems. The statistical approach to value 
statistics is meaningful since two-point phase space correlations of WF vanish (for 
distinct points) in the semi-classical limit as shown in a longer paper. We have 
proven that random states corresponding to the regime of quantum chaos have Gaus- 
sian WF value statistics, with the relative standard deviation divided by the mean 
value diverging as a/W ~ s/N in the semi-classical limit N ^ oo. In addition we 
have shown that WF statistics in a classically fully chaotic system, initially in coher- 
ent state, tends to a Gaussian statistics within the Ehrenfest time scale tr- We believe 
that our results may shed some light onto the related studies of decoherence^^' 
and the fidelity of two initially equivalent states which undergo two slightly different 
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Fig. 4. Percentage of phase space with negative WF as a function of scaled time for the quantized 
saw-tooth map, K = 0.5, at several different dimensions N. 



time-evolutions.^) The randomness of WFs in classically chaotic systems suggests 
a description by means random functions. The latter is the subject of the more 
detailed forthcoming publication 13). 
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